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Seasonal  water-temperature  changes  influence  stage-discharge  and  depth- 
discharge  /relations  in  rivers  and  in  conveyance  channels  and  naviga- 
tion channels  where  wateij/discharge  is  constant.  The  major  influent 
/is  related  to  changes  in  bed  configuration  which  accompanies  water-/ 
temperature  changes. 

Stability  analysis  provides  a theoretical  basis  for  predicting  bed- 
configuration  changes,  and  using  the  theoretical  relations^with  one 
empirical  coefficient  defined  from  field  data,  graphical  relations  are 
developed  to  use  as  design  curves  for  estimating  the  velocity  and  flow 
depth,  given  the  design  unit  water  discharge,  bed  material  size  distri- 
bution, water  temperature,  and  channel  slope.  The  design  curves  are  for 
channels  with  bed  materi^s  from  0.1  to  0.3  mm  median  diameters,  water 
temperatures  for  0 to  30^C,  and  slopes  from  0.00005  to  0.0012. 
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INTRODUCTION 


Field  observations  show  that  for  some  natural  streams  and  con- 
veyance channels  flowing  at  constant  discharge,  a change  in  water 
temperature  alone  will  promote  a drastic  change  in  bed  configuration 
(Colby  and  Scott,  1965;  U.S.  Army  Corps  of  Engineers,  1969;  Culbertson 
and  others,  1972).  Until  recently,  there  was  no  theoretical  basis  to 
identify  the  mechanism  responsible  for  the  changes  in  bed  configurations 
or  to  predict  the  effects  of  temperature  upon  bed  forms  and  flow  resist- 
ance. Despite  intense  study  for  many  decades  the  basic  hydraulic  and 
transport  processes  in  alluvial  channels  can  be  described  only  in  an 
approximate  manner. 

A major  advance  in  the  theory  of  alluvial  channel  flows  was  the 
application  of  linear  stability  analysis  to  identify  the  regions  of  flow 
for  which  the  stream  bed  would  be  stable;  that  is,  the  flow  conditions 
under  which  a small  disturbance  of  an  initially  plane  bed  would  be 
dampened.  Early  theories  of  Exner  (1925)  and  Polya  (1937)  predicted 
instability  of  the  bed,  but  no  stable  regions.  Anderson  (1953)  first 
pointed  out  the  importance  of  the  strong  interaction  between  water 
surface  and  bed  waves  for  antidunes,  bed  and  water  surface  configura- 
tions that  had  been  documented  many  years  earlier  by  Pierce  (1916)  and 
Gilbert  (1914).  However,  it  was  Kennedy  (1961,  1963)  who  first  put 
forth  a comprehensive  theory  to  account  for  dunes,  a plane  bed,  and 
antidunes  by  means  of  linear  stability  analysis.  Kennedy  used  a poten- 
tial flow  model  with  an  empirical  relation  between  flow  velocity  and 
sediment  transport;  in  later  studies  the  basic  theory  was  extended  to 
the  formation  of  ripples  in  closed  conduits  (Kennedy,  1964)  and  to 
ripples  formed  in  oscillatory  flows  (Kennedy  and  Falcon,  1965).  An 
excellent  review  article  summarizes  the  basic  theory  and  its  appli- 
cation (Kennedy,  1969)  . 

Many  others  have  contributed  to  the  stability  approach  (Reynolds, 
1965;  Gradowczyk,  1968;  Smith,  1970;  Engelund  and  Hansen,  1966;  Hayashi, 
1970;  Engelund,  1970),  and  it  is  well  established  that  stability  analy- 
sis can  explain,  at  least  to  a first,  approximation,  the  essential 
features  of  the  observed  data  on  bed  configurations  in  alluvial  chan- 
nels. 


The  earlier  stability  analyses  of  Kennedy  (1963)  and  Reynolds 
(1965)  were  based  on  potential  flow  models  and  so  were  not  capable  of 
including  any  dynamic  effects.  Engelund  (1970)  introduced  a real  fluid 
model  that  provided  a more  realistic  description  of  the  dynamic  forces, 
the  flow,  and  the  coupled  sediment  transport.  However,  all  of  the 
linear  stability  analyses  have  in  common  two  major  shortcomings. 

First,  because  of  the  linearity,  only  first  order  effects  are  con- 
sidered, and  this  can  serve  to  identify  only  the  initial  instabilities. 
The  equilibrium  height  and  final  wavelength  to  which  dunes  can  grow 
cannot  be  predicted  by  linear  methods,  nor  can  the  flow  separation  or 
nonsymmetry  of  the  bed  forms  be  described.  Second,  although  the  stabil- 
ity analyses  have  provided  a great  deal  of  insight  into  the  mechanism  of 


flow  in  alluvial  channels,  they  have  not  led  to  any  practical  advances 
in  design  of  stable  channels  or  in  prediction  of  bed  forms  and  flow 
resistance.  If  one  knows  all  the  flow  and  sediment  variables  required 
for  the  stability  analyses,  all  the  essential  features  of  the  flow 
required  for  a practical  design  problem  are  known. 

Recently,  Engelund  and  Freds^e  ( 1 9 7 A ) introduced  some  important  new 
concepts  into  stability  analysis.  They  began  with  the  vorticity  trans- 
port equation  to  describe  the  flow,  following  the  basic  model  proposed 
by  Engelund  (1970),  and  they  proposed  that  the  stability  of  the  bed, 
described  by  the  imaginary  part  of  the  bed-wave  migration  velocity,  was 
dependent  upon  the  ratio  of  suspended  load  to  bed  load.  A temperature- 
dependent  transport  relation  in  which  it  is  assumed  that  all  particles 
smaller  than  a critical  size  will  enter  into  suspension  was  introduced, 
and  because  the  critical  size  depends  on  the  ratio  of  particle-fall 
velocity  to  shear  velocity,  the  analysis  predicts  very  well  the  tran- 
sition from  dune  to  plane  bed  with  decreasing  water  temperature. 

The  purpose  of  this  report  is  to  extend  the  analysis  presented  by 
Engelund  and  Freds^e  (1974)  to  a greater  range  of  conditions  so  that 
stability  boundaries  between  dunes  and  a plane  bed  can  be  identified  and 
presented  graphically  as  aids  in  practical  channel-design  problems.  In 
the  following  section,  some  of  the  basic  elements  of  stability  analysis 
as  developed  by  Engelund  and  Freds^e  (1974)  are  presented.  This  is 
followed  by  a summary  of  the  results  of  the  stability  analysis,  pre- 
sented graphically,  and  finally,  it  is  shown  how  an  extension  of  the 
basic  equations  can  provide  an  estimate  of  mean  velocity  and  flow  depth 
if  the  channel  width,  design  discharge,  particles  size  and  water  tem- 
perature are  known.  A numerical  example  is  given  in  the  appendix. 


STABILITY  ANALYSIS 


1.  General  Approach 

The  usual  approach  to  fluid-stability  analysis  is  to  introduce  a 
small  perturbation  to  the  flow  and  to  investigate  how  the  amplitude  of 
the  perturbation  increases  or  decreases  with  time.  Usually,  the  per- 
turbation takes  the  form  of  a sinusoid  so  that  its  amplitude  can  be 
described  in  complex  form.  If  the  imaginary  part  grows  with  times,  the 
flow  is  unstable,  while  if  the  imaginary  part  is  dampened,  the  flow  is 
stable. 

For  flow  over  a sand  bed,  the  stability  analysis  is  complicated 
because  the  flow  equation  and  the  sediment  transport  must  be  coupled  and 
because  the  model  must  account  both  for  the  dune  case  where  the  water- 
surface  and  bed  waves  are  out  of  phase  and  for  the  antidune  case  when 
the  bed  and  water-surface  waves  are  in  phase.  The  migration  of  the  bed 
forms  must  be  included  and  transport  both  of  bed  load  and  suspended  load 
has  to  be  considered. 
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In  the  following  section,  the  stability  theory  derived  by  Engelund 
and  Freds^e  (1974),  that  accounts  for  all  these  factors  Is  outlined. 

The  discussion  here  follows  very  closely  the  original  work  of  Engelund 
and  Freds^e  (1974),  where  a complete  derivation  of  the  stability  theory 
is  given. 

2.  Summary  of  the  Theory 
(1)  Flow  equations 

The  perturbed  flow  for  two  dimensional  flow  over  a slightly 
sinusoidal  bed  may  be  described  by  the  vorticity  transport  equation 
(Engelund,  1970) 


dm 

di  = 


eV2u > 


(1) 


in  which  e is  the  eddy  viscosity.  The  vorticity  w is  defined  by 


3 u 

3 y 


(2) 


where  x is  a coordinate  in  the  mean  flow  direction,  y is  a vertical 
coordinate,  and  u and  v are  the  velocity  components  in  the  x and  y 
directions,  respectively.  The  stream  function  4*  for  the  flow  is  defined 
by 


u = 


3^ 

3 y 


and 


(3) 


If  Eqs.  (2)  and  (3)  are  substituted  into  the  vorticity  transport 
equation,  we  get 


3»  3(72»;  3fV2ifrj  3(72a»;  _ (4) 

3 y 3x  dx  3 y H e ^ 


3 


Next,  the  distribution  of  the  suspended  load  in  a nonunlform  flow 
Is  described  by  a local  continuity  equation  expressing  the  equilibrium 
between  settling,  diffusion  and  convection  of  sediment 


dC 
It  ' 


W + eV2C 

*y 


(5) 


where  C is  the  concentration  of  the  water-sediment  mixture  and  ^ is  the 
fall  velocity  of  a representative  sediment  particle.  When  the  uniform 
flow  is  perturbed,  it  is  assumed  that  the  perturbation  of  the  flow  is 
periodic  and  is  described  by  the  following  relations  (Engelund  and 
Freds^e,  1974), 


h = hoexp  j ikD  (§  - | t)J 

(6) 

n = r\Qexp  j ikD  - | i)j 

(7) 

a = i 

(8) 

<|»  - 'i>0  + 'P 

(9) 

* = {ikD  (f- ft)) 

(10) 

c = c + c 

(ID 

o 

c = ~ f2(^)  exp  {iWJ(|  - | t)| 

(12) 

In  which  h is  the  local  height  of  the  sand-bed  wave  (see  definition 
sketch.  Fig.  1),  hQ  is  the  amplitude  of  the  sand-bed  wave,  kD  is  the 

dimensionless  wave  number  with  k = 2-n IL,  L is  the  wave  length,  a is  the 
complex  migration  velocity  of  the  bed  waves,  n is  the  local^height  of 
water-surface  wave,  4>0  is  the  undisturbed  stream  function,  ^ is  the 
disturbed  3tream  function,  Uj^0  is  the  finite  wall  slip-velocity,  CQ 
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and  C are  the  undisturbed  and  the  disturbed  sediment  concentrations, 
and  /x(§)  and  f2ty)  are  unknown  functions  to  be  found  by  solving  the 
vorticity  transport  equation  (4)  and  sediment  continuity  equation  (5). 


Figure  1. — Definition  sketch. 


To  solve  Eqs.  (4)  and  (5),  the  undisturbed  velocity  and  concentra- 
tion distributions  are  required.  The  undisturbed  velocity  profile  can 
be  described  by  the  parabolic  velocity  distribution  in  the  major  part  of 
the  flow  (outside  the  viscous  sublayer)  from  integration  of  the  flow 
equation.  Assuming  a constant  value  of  the  eddy  viscosity  e,  the  un- 
disturbed velocity  profile  is  described  by 


U = U 


so 


(13) 


where  U and  U ± are  the  undisturbed  surface  velocity  and  shear 

so  *o 

velocity,  respectively,  and  e is  given  by  (Hinze,  1959) 


e = 0.077  U*  D (14) 

*o 

The  quantity  U A is  actually  the  effective  shear  velocity  V\  over  a 

plane  bed  and  is  found  from  a graphical  solution  of  the  Einstein- 
Barbarossa  equation  (Fig.  2,  from  Vanoni  and  Brooks,  1957).  In  order 
to  obtain  a simple  analytical  solution,  the  equation  of  a parabola  (13) 
is  replaced  by  a cosine  given  by  the  expression 

U = Uso  COS  [BC1  - £)]  (15) 
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u * 

with  [V  = 14  jj — to  approximate  the  parabola.  A finite  wall  slip 

iso 

veloc:  Lty  must  be  introduced  with  the  assumption  of  constant  eddy 
viscosity  (Er.gelund,  1970) 


U,  = £/  cos  8 
bo  so 


(16) 


which  is  found  by  matching  the  velocity  profile  outside  the  viscous 
sublayer  with  a known  logarithmic  velocity  profile  near  the  bottom. 
For  a hydraulically  rough  bed,  the  matching  condition  is 


Uh 

jf2-  = 1.9  + 2.45  In  j—  = K 
*o  s 


(17) 


where  if  is  a constant  and  k is  the  representative  sand  roughness  and  is 


taken  to  equal  2.5  with  d ^ the  mean  diameter  of  the  bed  material. 


The  undisturbed  concentration  can  be  found  by  solving  Eq.  (5)  in 
the  steady  uniform  flow  as 


Co  = Cbo  ?) 


(18) 


where  is  a sediment  concentration  at  the  bed  level. 

When  Eqs.  (9)  through  (12)  are  substituted  into  Eqs.  (4)  and  (5), 
we  obtain  the  following  pair  of  ordinary  differential  equations  in  / 

and  on  the  assumptions  that  h « D and  a « U: 


n-  + 


e 

JWU 


■pint 


(19) 


and 


n + ft'i  - [lkD>z  * ikD  rl  h - m ar 

l j y 


(20) 


where  £ = y/p  and  all  differentiations  are  made  with  respect  to  £ . 

y y 
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Analytical  solutions  to  Eqs.  (19)  and  (20)  are 


The  coefficients  6'2,  Cy  C , and  a 2 in  Eqs.  (21)  and  (22)  can  be 

found  from  six  boundary  conditions  described  in  Sec.  2—  (2 ) . Once  the 
functions  and  /2  are  determined,  the  migration  velocity  of  bed  waves 

is  obtained  Dy  solving  the  sediment  continuity  equation 


9x 


Cl  - m) 


Zh 

U 


or  in  another  form 


a 


11  - m!  ^ - 


(23) 


(24) 


where  m is  the  porosity  of  the  sand  bed  and  q^,  is  the  total  bed  material 
load  which  is  the  sum  of  the  bed  load  q ^ and  the  suspended  load  qg.  The 
value  of  the  imaginary  part  of  the  migration  velocity  determines  the 
stability  of  the  sand  bed.  Eq.  (6)  can  be  written  as 


h = hQ  exp  (had;)  exp 


(25) 


where  ka  ■ is  called  the  amplification  factor,  if  a.  is  found  to  be 

positive,  the  sand  bed  will  be  unstable  because  the  amplitude  of  the 
perturbation  will  increase  in  time.  Correspondingly,  the  bed  will  be 
stable  if  < 0.  Furthermore,  the  real  part  a determines  the  migra- 
tion velocity  of  the  bed  wave  in  the  flow  direction  as  well  as  the  type 
of  the  bed  form.  If  > 0 , the  bed  wave  migrates  downstream  and  the 

bed  is  about  180°  out  of  phase  with  the  water  surface.  The  sand  bed  has 
the  characteristics  of  dunes.  In  the  case  of  antidunes,  bed  and  water 
surface  are  approximately  in  phase  and  < 0 . 

(2)  Boundary  conditions 

To  obtain  the  four  coefficients  £7  , (?2,  Cy  and  C ^ in  Eq.  (21), 
the  following  four  conditions  are  employed: 

(i)  The  kinematical  condition  at  the  sand  bed  is 


dh  _ dh 
dt  9 1 


9 h 

dX 


= V 


(26) 
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at  r = h/l).  When  eqs.  (3),  (6),  (9)  and  (10)  are  Introduced,  Eq.  (26) 

7/ 

Is  reduced  to 

f^O)  = 1 (27) 

assuming  a/7/^  « 1 and  h/D  « 1. 

(ii)  The  same  kinematical  condition  at  the  surface  yields 


U T) 

f (x)  = _«£  _£ 

Jl(i;  u,  h 

bo  o 


(28) 


and  the  shear  stress 


T ,3m  , 3u  , „ 2 

— = e (r—  + — ) = US 
p dy  dx  * 


(29) 


must  vanish  at  the  free  surface.  Eqs.  (28)  and  (29)  lead  to  the  second 
boundary  condition 


f'{(  1)  + fx(  1) 


(30) 


(iii)  The  normal  stress  must  vanish  at  the  free  surface.  To  a 
first  order  approximation,  this  is  expressed  as 


P 

s 

P 


= 2 


(31) 


The  surface  tension  is  neglected  and  P is  the  water  pressure 

s 

just  below  the  surface.  To  eliminate  P from  Eq.  (31),  the  equation 

s 

of  motion  at  5 = 1 is  considered 


du 

di 


3_ 

3x 


* f) 


+ eV2  u + gS 


(32) 
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in  which  S is  the  mean  slope  of  the  channel  equal  to  the  slope  of 
the  undisturbed  water  surface.  must  satisfy  the  relation 


e 


+ g S 
J o 


= 0 


(33) 


Based  on  Eqs.  (31),  (32),  and  (33),  the  boundary  condition  becomes 

r u d 1 

f{"( v + I -ikD  ~~  - 3(kDJ2\  f^fD  + gL  itcDf  (1)  = o (34) 

80 


(iv)  The  condition  (17)  is  also  assumed  to  be  valid  in  the  per- 
turbed flow,  i.e. 


U. 

b _ V 

u*  ~K 


(34a) 


at  C7  ~ h/D.  Determining  the  local  value  of  U A by  Eq.  (29),  the  linear 

y 

form  of  Eq.  (34)  is 


2U.D  Ml  D 

n<°>  - -wr  nl0>  * <kB>1  h(0>  - f - tpI - -i r1  * 0 (35) 

bo 


with  all  differentiations  being  with  respect  to  K . 

y 

Two  boundary  conditions  are  required  to  obtain  the  two  constants 
a ^ and  a ^ in  Eq.  (22). 

(k)  The  relation  between  the  bed  concentration  C,  and  the  local 

b 

value  of  the  shear  stress  is  assumed  to  follow 


V2 


(36) 


where  iq  is  the  undisturbed  bed  shear  stress.  After  the  substitution, 
the  linearized  form  turns  out  to  be 
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C + h 


dC  1C, 

o _ h 

du  U 


, 3Y  . 2 dll  . 

f — T — + rl  7 — ) 

dy 


bo 


or 


f2(0)  = 3cbo  [If  ■ fi(0)  + 8 Un  B] 


(37) 


(ii)  At  the  water  surface,  the  vertical  sediment  flux  must 
vanish 

c ~ + WC  = 0 


Substitution  of  Eq.  (12)  for  C into  the  above  equation  yields 


ftfl)  + ~ f2(l)  = 0 


(38) 


(3)  Sediment  transport 

To  determine  the  migration  velocity  of  the  bed  wave  from  Eq.  (24), 
knowledge  of  sediment  transport  is  needed.  Sediment  transport  rates 
are  affected  by  a large  number  of  factors,  and  there  is  still  no  method 
available  at  present  for  accurately  estimating  the  sediment  discharge 
under  different  flow  conditions.  In  this  analysis,  the  bed  load  trans- 
port rate  q ^ is  assumed  to  follow  the  Meyer-Peter  formula 

% 3/ 

= - = 8 ( 9-0.047 ) /2  n<n 

b /(e-l)g  db 3 


where  s is  the  relative  density  of  the  sediment,  g is  the  accelera- 
tion of  gravity,  d ^ is  the  mean  diameter  of  the  grain  moving  as  bed 

load,  and  9 is  a parameter  representing  the  ratio  between  the  effective 
shear  stress  at  the  bed  and  the  stabilizing  forces  on  a sediment  par- 
ticle in  the  bed 


9 


p(s-l)g  db 


+ \iS 


(40) 
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where  x'  is  the  effective  shear  stress,  S is  the  bed  slope  and  u is  a 
dynamic  friction  coefficient  estimated  to  be  0.1.  The  ratio  between 
disturbed  and  undisturbed  bed  load  transport  rates  is  then 


% = (9-0.047) 

%o  (e  -0.047; 

O 


(41) 


where 


0 = 


+ ViS 


o (s-l)g  d^  r~o 


in  which  Sq  is  the  undisturbed  plane  bed  slope. 

The  rate  of  suspended  sediment  transport  q is  found  from 

O 


' D+  n 


qs  = 


uCdy  » qso  - UhoCbh  + 


where 


s:  k;  * 


V hCo  * °V2V 


b 


(42) 


1 so 


= D\ 

J o 


VC  dl 
o ''y 


is  the  suspended  load  transport  rate  over  an  undisturbed  flat  bed. 

The  ouantities  q^  and  qgQ  can  be  estimated  from  Eqs.  (39)  and 
(42)  provided  d ^ and  C are  known.  These  two  quantities  can  be 
evaluated  based  on  the  following  two  assumptions: 

(i)  Sediment  particles  finer  than  a certain  critical  size  all 
enter  into  suspension  while  larger  particles  move  as  bed  load.  A 
preliminary  criterion  for  determining  the  critical  particle  size  was 
suggested  by  Engelund  (1973) 


w = 0.8  ur 

a * 


(43) 
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Once  the  critical  settling  velocity^  and  the  corresponding  size 

are  determined,  the  median  size  of  the  bed  load  can  be  estimated  from 
the  bed  material  size  distribution. 


(ii)  No  segregation  takes  place  so  the  moving  particles  keep  the 
original  composition  of  the  bed  material.  This  assumption  states  that 


b _ (total  volume  of  particles  transported  as  bed  load) 

V ~ (total  volume  of  particles  transported  as  suspended  load) 

s 

Percent  bed  material  coarser  than  d 


Percent  bed  material  finer  than  d 


(44) 


where  in  a unit  area 


Vb  = °'H  • 


and 


S 


V = I c dy  - DC. 


8 >d  ° 


bo 


[¥]'. 


W ^ is  the  fall  velocity  of  the  representative  size  of  suspended  load 

d and  percent  loads  are  obtained  from  assumption  (i).  Eq.  (44)  yields 

s 

the  value  of  C,  . 

bo 

It  is  interesting  to  note  the  temperature  effects  on  the  sediment 
transport  rates.  The  primary  effect  of  temperature  is  its  large  in- 
fluence on  the  viscosity  of  the  fluid  and  on  the  fall  velocity  of  sand 
particles  in  suspension.  A decrease  in  the  water  temperature  increases 
the  fluid  viscosity  and  decreases  the  fall  velocity  of  suspended  par- 
ticles. Based  on  assumption  (i),  the  decrease  in  temperature  increases 
the  critical  size  and  thus  increases  the  percent  suspended  load  of  the 
total  load.  Then  from  assumption  (ii),  this  will  increase  the  value  of 

C,  and  will  cause  the  concentration  distribution  to  be  more  uniform. 

bo 

An  example  of  temperature  effects  on  the  sediment  transport  rates  cal- 
culated from  Eqs.  (39)  and  (42)  is  shown  in  Fig.  3.  The  decrease  in 
temperature  significantly  increases  the  suspended  load  discharge  but 
has  only  slight  effect  on  the  bed  load.  In  this  example,  the  velocity 
and  the  depth  are  assumed  unchanged  with  temperature.  However,  if  the 
flow  conditions  are  close  to  the  transition  zone,  the  temperature 
changes  may  change  the  bed  form  and  its  resistance  to  flow.  In  this 
case,  the  decrease  in  temperature  will  increase  the  velocity  and 
decrease  the  depth  for  a given  flow  discharge,  and  therefore  tremen- 
dously increase  the  ability  of  the  stream  to  carry  the  sediment  load. 
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Figure  3. — Temperature  effects  on  the  sediment  transport  rates 
(q  = 50  cfs/ft,  = 0.21  mm,  a = 1.23,  S = 0.00015  and  F * 0.25) 

(4)  Migration  Velocity 


The  migration  velocity  determines  the  stability  of  the  bed  and  also 
the  form  of  bed  roughness  as  stated  in  Sec.  2-(l).  Its  value  can  be 
determined  from  the  sediment  continuity  equation  (24).  If  the  part  of 
the  total  sediment  discharge  carried  in  suspension  is  converted  to  an 
equivilant  bed  migration  velocity;  Eq.  (24)  can  be  written  as 


(ab 


V 


(1-m) 


dh 

3x 


dJb  + *_% 

3X  dX 


(45) 


where  a ^ and  are  components  of  migration  velocity  affected  by  bed 

load  and  suspended  load  transport  rates,  respectively.  Eq.  (45)  can  be 
written  correspondingly  as 


a,  (1-m)  ~ 
b 3a; 


3X 


and 


a (1-m) 
8 


3 h 
3x 


3x 


(46) 


(47) 


15 


Substituting  Eqs.  (41)  and  (42)  fojr  a,  and  q into  Eqs.  (46)  and  (47) 

D ~ 8 

yields  the  linear  form 


ab  1 -m  (0  -0.047)  (3-l)g  duD 

o u b 


1 DO 


gd,(8- 1) 

-f^(0)  +6  tan  6 - ikD  y 


*1 


and 


*.  - rs  {£|>W  + " W]  Vi.) 


(48) 


(49) 


When  values  of  a,  and  a are  calculated  from  Eqs.  (48)  and  (49),  the 

b s 

equivalent  total  migration  velocity  is 


a = a,  + a 
b s 


(50) 


The  imaginary  part  becomes 


ai  = abi  + asi 

t-  U U O Is 


(51) 


The  values  of  as  strongly  depend  on  the  ratio  between  bed  load  and 
suspended  load.  In  Fig.  4,  an  example  of  the  quantities  a ^ and 

is  presented.  The  figure  shows  that  a . < 0 and  a,  . > 0 when  the 

st  bx 

Fronde  number  is  less  than  0.8.  This  indicates  that  in  the  lower  flow 
regime  the  suspended  load  stabilizes  the  bed  while  the  bed  load 

reverses  the  effect.  Both  a,  . and  -a  . increase  with  an  increase  in 

bt  si 

Froude  number.  However,  when  the  Froude  number  is  larger  than  a cer- 
tain value,  -a  . increases  much  faster  than  a,  . . Hence,  when  the 

8t  bt 

Froude  number  is  low,  a • is  positive  indicating  bed  instability  with 

t 

bed  forms  of  dunes,  while  for  high  Froude  number,  a ^ becomes  negative 
indicating  bed  stability  (flat  beds). 

In  Fig.  4,  the  temperature  effects  on  a,  . and  a . also  are  shown. 

A decrease  in  temperature  increases  both  a , . and  -a  . • Nevertheless, 

bt  st 
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Figure  4. — Temperature  effects  on  the  imaginery  migration  velocity  of 
bed  wave  (D  = 0.33  ft,  = 0.21  mm,  a = 1.23,  kp  = .21) 

the  increase  in  -a  . is  larger  because  of  the  large  temperature  effect 

O t' 

on  the  suspended  load  as  shown  in  Fig.  3.  This  indicates  that  dunes 
disappear  at  lower  Froude  numbers  when  the  temperature  is  low. 

3.  Results  of  the  Stability  Analysis 

As  shown  in  Fig.  4,  for  a given  wave  number  kD , depth  D,  median 
size  of  bed  material  d gradation  coefficient  of  sediment  size  dis- 
tribution a and  water  temperature  T,  the  bed  may  be  stable  or  un- 
stable depending  on  the  flow  velocity.  When  kD  is  plotted  against 
Froude  number  F keeping  other  variables  unchanged,  a set  of  transition 
curves  distinguishing  dunes  and  plane  bed  is  obtained  as  shown  in 
Fig.  5.  In  the  figure,  the  instability  of  the  lower  flow  regime  dis- 
appears at  a Froude  number  about  0.84  and  temperature  30°  Celsius.  As 
pointed  out  in  Sec.  2- (4),  this  transition  between  dunes  and  plane 
bed  depends  on  the  ratio  between  bed  load  and  suspended  load.  Because 
temperature  has  significant  effects  on  this  ratio,  the  decrease  in  water 
temperature  lowers  the  stability  limits  of  the  transition  line  and  thus 
may  remove  dunes  from  the  sand  bed.  From  the  figure,  it  is  concluded 
that  the  transition  from  dunes  to  plane  bed  always  takes  place  at  very 
large  wavelengths. 


I 


Figure  6 shows  more  specifically  the  temperature  effects  on  bed 
stability.  The  curves  represent  the  upper  limit  for  the  existence  of 
dunes  as  functions  of  D,  T,  a and  mean  flow  velocity  V.  For 

a given  depth  and  velocity  if  the  temperature  decreases  to  a certain 
value  the  dunes  disappear.  For  example,  when  D = 10  ft  and  V = 6.5 
fps,  if  T Is  less  than  30°  Celsius,  the  bed  form  changes  from  dunes 
toward  a plane  bed.  At  a still  lower  temperature,  the  change  from 
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Figure  5 . r-Tenroerature  effects  on  the  instability  regime 
* (D  = 0.33  ft,  d50  = 0.3  mm,  o = 1.23; 

dunes  to  plane  bed  occurs  at  a lower  velocity.  However,  if  the  flow 
conditions  are  far  below  the  transition  zone,  the  temperature  effect  is 
not  significant. 

The  most  obvious  example  of  temperature  effect  is  that  reported  for 
the  Missouri  River  (U.S.  Army  Corps  of  Engineers,  1969),  where  a de- 
crease in  temperature  during  the  autumn  first  causes  a stretching  of  the 
dunes  and  finally  smoothes  out  the  bed  so  that  it  becomes  completely 
plane.  When  Missouri  River  data  (see  Table  1)  are  plotted  on  Fig.  6, 
they  fall  mostly  between  the  transition  curves  of  T = 20°  and  10° 
Celsius.  Then,  as  indicated  in  this  figure,  the  decrease  in  temperature 
lowers  the  critical  curve,  changes  the  bed  form  from  dunes  to  plane  bed, 
decreases  the  depth  and  increases  the  velocity. 


f 

li 


FLOW  RESISTANCE 


When  the  flow  parameters  nearly  satisfy  the  conditions  for  transi- 
tion between  dunes  and  plane  bed,  the  effects  of  temperature  on  flow 
characteristics  become  significant.  The  decrease  in  temperature 
smoothes  out  the  sand  bed  and  therefore  reduces  its  resistance  to  flow. 
This  reduction  in  flow  resistance  increases  the  flow  velocity  and  de- 
creases the  depth  for  a given  unit-width  flow  discharge  as  shown  in 
Fig.  6.  Nevertheless,  the  amount  of  change  cannot  be  evaluated  from 
Fig.  6.  An  attempt  to  predict  the  temperature  effects  on  resistance  to 
flow  is  initiated  here  based  on  stability  theory. 


1.  Theoretical  Analysis 


The  shear  stress  can  be  related  to  the  velocity  gradient  as 


T ,WU  , dV . 
p = + 


(52) 
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MEAN  VELOCITY,  V,  IN  FEET  PER  SECOND 


A 


1 1 i 


In  practice,  the  bed  shear  stress  is  commonly  composed  of  the  surface 
drag  and  the  form  drag. 


T = T ' + f' 


(53) 


Similarly,  b>  applying  stability  theory,  the  shear  stress  is  composed 
of  undisturbed  and  disturbed  shear  stresses. 


T = + T 


(54) 


where  tq  is  the  undisturbed  bed  shear  stress  acting  on  the  flat  bed 

under  a uniform  flow  field  and  is  due  to  the  grain  roughness  of  the 
bed.  This  indicates  that 


To  = T 


(55) 


and  that  the  form  drag  is  mainly  caused  by  the  disturbed  flow  field 
due  to  the  bed  forms. 

To  eliminate  the  effects  of  local  nonuniformity,  the  bed  shear 
stress  is  evaluated  from  averaged  flow  characteristics  over  a wave 
length  as 


T = T + X 

O 


where  L is  the  bed  wave  length,  and  1 is  a correction  factor  to 
relate  the  energy  dissipation  due  to  form  roughness  to  the  disturbed 
flow  field  on  the  bed.  The  disturbed  shear  stress  is  averaged  over  a 
quarter  of  the  wave  length  and  then  multiplied  by  2 to  account  for  the 
flow  disturbance  caused  by  the  bed  wave  of  amplitude  (2h^) . Comparing 
Eqs.  (53)  and  (54)  leads  to 


r „ rU 4 

> 

o' ' 

1 

(56) 

- - > (i  f 4 T *) 


(57) 


\pplying  the  stability  theory,  the  quantity  t can  oe  expressed  by 


Substituting  F.q . (10)  for  l1  Into  Eq . (58),  the  disturbed  shear  stress 
In  steady  nonuniform  flow  Is 


~ = -2 U\o  ~~  | f'(0)  - 8 tan  Bj  exp  <i2v  (59) 

This  equation  can  be  rewritten  with  the  aid  of  Eq.  (48)  as 


t 1 \a-m)  J(s-Y)g  db  (* 

p 12  db  f reo-0.047;  \ ab 


+ ^(eQ~0. 047  J (S-I)g  d^exp  fillip)  (60) 

Substituting  Eq.  (60)  into  Eq.  (57)  for  t and  then  integrating  Eq. 

(57),  the  real  part  of  the  resulting  relation  equals  the  form  drag, 
that  is 


si  * J<s-1>3  db  ,2ho>  r 

p ” 6ir  yj  (0^0.047)  dh  l br 


(61) 


where  without  loss  of  generality  m is  assumed  to  be  zero.  The  tem- 
perature effect  on  the  form  drag  is  clearly  shown  in  Eq.  (61).  The 
decrease  in  temperature  stretches  the  dunes  and  reduces  the  migration 
velocity  and  the  amplitude  of  the  bed  forms.  Then  from  Eq.  (61), 
these  effects  decrease  the  form  roughness  of  the  bed. 

To  check  the  validity  of  Eq.  (61),  field  data  from  the  Missouri 
River  and  the  Rio  Grande  (see  Table  1)  are  employed  to  determine  the 
correction  factor  X.  The  quantities  t",  hQ,  kD , D and  0^  are  obtained 

from  the  field  data,  and  the  values  of  abr>  a£>£»  an<*  are  evaluated 

from  the  stability  theory.  The  results  are  plotted  in  Fig.  7 to  relate 
X and  T.  It  is  found  that  X decreases  with  temperature  and  approaches 
unity  at  low  temperatures.  The  reasons  are: 

(i)  At  high  temperatures,  the  bed  configuration  is  dunes  that 
have  approximately  triangular  profiles  with  long  gentle  upstream 
slopes  and  short  steep  downstream  slopes  causing  large  separation  zones 


- a 


bi 


1.2 (kD) 
D 


J( eo-0. 047 )7e^l)g  d*  ] 


r 

[! 


In  the  flow.  These  profiles  and  large  separation  zones  of  dunes  are 
not  well  simulated  by  the  stability  theory,  which  assumes  a sinuous 
bed  profile  and  small  wave  amplitude  and  separation  zone.  Therefore, 
the  energy  dissipation  caused  by  dunes  is  expected  to  be  larger  than 
estimated  from  the  theory.  The  correction  factor  accounts  for  this 
difference. 

(ii)  At  low  temperatures,  the  bed  wave  is  smoothed  out  approach- 
ing a sinuous  profile  as  described  by  Eq.  (6).  Thi3  enables  the  theory 
to  simulate  the  flow  field. 


Figure  7. — Temperature  effects  on  the  resistance  correction  factor. 

The  scatter  in  Fig.  7 is  mainly  caused  by  fluctuations  in  the 
amplitudes  of  the  bed  waves.  By  using  a more  consistent  method  to 
analyze  the  bed  forms,  the  scatter  could  be  reduced. 

2.  Design  Principles 

Eq.  (61)  involves  properties  of  the  bed  forms  and  is  not  practical 
for  design  purposes.  A more  generalized  method  is  thus  developed. 

Eq.  (61)  is  rearranged  in  the  form 


T (db  V)  J<‘-Vg  *50  (,  ho  \abf 

1‘U1' 


*bi 

v 


L. 


1 . 2 (kUl 

VD 


jreo-0.0t7)(e-i)ff  d*  Jf  = {/(yip: 


-i )g  d 


32 


-0.047 


° v (62) 


where  9^  = (s-\)g  d — + and  S ^ is  the  mean  bed  slope.  The 

coefficient  o is  related  to  the  disturbed  flow  characteristics  and 
/ is  a function  of  size  distribution  and  flow  conditions.  In  Eq.  (62), 
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Figure  8b. — Relation  for  disturbance  coefficient  (d5o  = 0.2  mm 

S = 0.000075) 
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Figure  8c. — Relation  for  disturbance  coefficient  (dsn  * 0.2  mm 


DISTURBANCE  COEFFICIENT 


DISTURBANCE  COEFFICIENT 


DISTURBANCE  COEFFICIENT 


10.0 


FROUDE  NUMBER 

Figure  8g. — Relation  for  disturbance  coefficient  (djg  = 0.2  mm, 

S = 0.0012) 
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the  value  \h  Id,  can  not  be  evaluated  from  stability  theory.  An 

o b 

alternate  method  is  used  for  determination  of  c as  follows: 

(i)  Given  q , F , T,  d^,  a,  and  S^,  compute  U \ from  Fig.  2. 

(ii)  Compute  based  on  the  Sec.  2- (3). 

(iii)  Compute  the  shear  velocity 


Uf  = gDSo  - (Uf) 


(iv)  Determine  the  quantity 


V 


(s-l)g  db 

0 -0.047 

o 


(v)  Assume  / = 0.1  and  d^  = 0.2  mm  and  compute  a from 


e = U'l2/( db  \ 

\ \eo-0.047  / 


The  calculated  results  of  c are  plotted  on  Fig.  8 as  functions  of 
q,  F,  T,  a and  S . These  figures  are  employed  to  determine  c for  the 

following  conditions: 

(i)  q = 5 - 300  (ft3/s)/ft. 


(ii)  SQ  = 0.00005  - 0.0012, 

(iii)  d^  = 0.2  mm. 


and  (iv)  T 


=0-30°  Celsius. 


The  values  of  decrease  with  temperature,  indicating  that  the 
decrease  in  water  temperature  will  reduce  the  degree  of  flow  disturb- 
ance by  bed  forms. 


To  determine  the  effect  of  particle  size  on  values  of  c,  families 
of  curves  defining  c as  a function  of  particle  size,  unit  water  dis- 
charge, and  water  temperature  can  be  constructed  for  constant  values  of 
slope  and  Froude  number.  A set  of  curves  for  slope  S = 0.00015  and 
Froude  number  F - 0.25  is  shown  in  Fig.  9.  The  relative  positions  of 
the  curves  are  not  changed  much  for  other  values  of  slope  and  Froude 
number,  so  the  curves  in  Fig.  9 can  be  used  to  adjust  for  bed  particle 
sizes  from  0.1  to  0.3  mm  the  values  of  a determined  from  Fig  8,  which 


? 1 !'■  'xarapi'.  of  the  adjustment  is 

given  in  tbt  tpp  ml  lx. 

• 1 : las.  ' mp  J i iea  1 coefficient.  Its 

value  a 1 'n.|  r ■ h by  determining  the  other 

quant i t i ■ In  Eq . (62)  for  given  field  flow  data  and  then  computing  /' 

!•  1 *•  ' '<■'  '’row  h’  • ’Mssis  -ippi  River,  the 

and  i.  to  Gra..;i . Conveyance  Channel,  which  show  sig- 
■ utill  d t develop  the  empirical 
re!  ■!  ' lh.  values  of  f are  tabulated  in  Table  1 and  are 

■I  ' • ’ unci  ions  oi  water  : emperatures  and  unit-width 

’I-  s of  ; i.ur. ‘ease  vith  decrease  in  temperature 
■ : cd  load.  figs.  2,  8,  9,  and  10  are 

• •. ip.  1 o flow  w!':h  churn*  teristius  close  to 

th  ' 1 it. r.  i • rnple  and  ;tej  by-step  procedure  are 

•<  • n.  * lo  • iv  line  flow  v Lo'  ity  and  depth  for  giver 

it  ’ _ir!  -ce  slope.  , temperature  T arid  median 

;i  '0  of  bed  materia)  i^r.> 


j 

! 


! 


An  G:>' 
Lure  evalu 

i l i ca 

III  i t -w  id ! 1 > 


•i  i . vi!  . -utter  in  the  i low  uh-ir  ir.t«ristlcs  with  tempera- 
! rves  is  in  Fig.  11.  it 

l.o  le  , :a  i tp  ture  fo  gi v<  i 

>e<  use  of  the  use  in  resist  aace  to  i l..>w 


- I ’ 
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1 
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v ’■ ■ fKWvTUiti;,  Ji.  °c 

• - Te-mri-s  •:  t 
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on  flo r real  stance 

•re/ft,  s 

1 002 , 

= 0 .2  m 0 - 3.3). 

1 1 . > 1 1 

Data 

fy  '.ho  pr 

ced  lieaign 

curves,  Figs.  2,  8, 

9,  and  10  are. 

t f p rod  no  • 

> .e  eid  da 

‘ 3 ’ mitej  in  Table  ) . 

T1  1 i ! 

' : h q it  i tie'  7,  C , T , and  d,--  The  calculated 

1 O JU  O.) 

I w . ">  and  the  velocity  V . The  calculated  depth 

i cl  v ill  : 1*  1 mured  vn luce-  as  shown  in  Fig.  12. 
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4.  Limitation  of  Design  Curves 

The  design  curves  for  evaluating  flow  resistance  have  several 
limitations: 


(i)  The  design  curves  are  developed  from  linear  stability  theory. 
The  basic  assumption  of  this  theory  is  that  hQ/D« 1.  Therefore,  for 

small  flow  depth,  say  D< 4 ft,  the  design  curves  are  not  applicable. 


(ii)  The  design  curves  are  only  applicable  to  flow  regimes  close 
to  the  transition  zone.  It  is  suggested  that  the  stream  power 

tV  = pgqS 

should  be  greater  than  about  0.27  (lbs/ ft) /s  to  apply  the  design  curves. 

(iii)  The  design  curves  cover  the  range  of  sediment  grain  sizes 
0.1  - 0.3  mm,  temperature*  0 - 30°  Celsius,  unit-width  flow  discharges 
5 - 300  (ft3/s)/ft,  and  slope  0.00005  - 0.0012.  Extension  of  the 
curves  outside  these  ranges  is  not  recommended. 
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Figure  12. — Comparison  between  measured  and  calculated  depths. 
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SUMMARY  AND  CONCLUSION 


When  the  flow  regime  is  close  to  the  transition  zone,  the  tem- 
perature effects  on  the  flow  characteristics,  sediment  transport  rate, 
bed  configuration  and  its  resistance  to  flow  become  quite  significant. 

It  is  found  in  many  cases  that  a decrease  in  temperature  removes  the 
dunes  from  the  sand  bed,  decreases  flow  resistance  and  depth,  and 
increases  flow  velocity  and  sediment  transport  rate. 

Stability  theory  that  includes  the  effect  of  temperature  on  the 
flow  characteristics  in  alluvial  channels  was  developed  by  Engelund 
and  Freds«5e  (1974).  The  theory  shows  qualitatively  the  temperature 
effects  on  the  flow  field.  To  quantitatively  predict  temperature 
effects,  the  idea  of  disturbed  flow  on  bed  forms  developed  from  the 
stability  theory  is  introduced  to  describe  the  form  roughness.  An 
equation  relating  the  form  drag  to  the  disturbed  shear  stress  is  thus 
derived.  This  relation  yields  suitable  form  drag  for  low  temperatures 
where  the  bed  is  smooth  but  not  for  high  temperatures  with  a dune  bed. 

A correction  factor  introduced  in  the  equation  indicates  that  the  form 
drag  is  about  twice  as  much  as  that  evaluated  from  the  stability  theory. 

A set  of  design  curves  is  developed  from  the  stability  theory 
to  estimate  the  flow  resistance  in  terms  of  evaluation  of  flow  depth 
and  velocity  including  temperature  as  a major  factor.  The  initial 
data  required  to  apply  the  procedure  are  the  unit  water  discharge, 
water  temperature,  energy  slope,  and  particle  size  distribution  of  bed 
material. 
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1. 


Numerical  Example 


A detailed  numerical  example  is  presented  to  show  the  application 
of  design  curves  Figs.  2,  8,  9,  and  10  to  determine  the  flow  resistance. 


The 

given 

information  is 

(i) 

<7 

= 

40  (ft3/s)/ft 

(ii) 

m 

i 

= 

5°C , 

(iii) 

s 

o 

= 

0.0002, 

(iv) 

dS  0 

= 

0.22  mm, 

and 

(v) 

d6S 

= 

0.25  mm, 

The  following  table  shows  the  calculation  procedure  with  explanations: 


F V U2  US  U1  'J\  D'  a f P U'f  D"  D q n 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10)  (11)  (12)  (13)  (14) 

0.275  4.59  1996  9.2xl08  28.5  0.161  4.03  0.197  .140  .247  .0313  4.86  8.89  40.80 
0.29  4.75  2066  1.02xl0928.5  0.166  4.28  0.164  .140  .239  .0260  4.04  8.32  39.52 


Col. 

(1): 

Col. 

(2): 

Col. 

(3): 

Col. 

(4): 

Co  1. 

(5): 

Col. 

(6): 

Col. 

(7): 

Col. 

(8): 

Assume  a Froude  number. 

^ 

Compute  the  mean  flow  velocity  from  V = F ' <gq)  ' . 
Compute  U2  = V/y/gd^S^. 

Compute  U3  = V^/gvS  . 

Find  U1  from  Fig.  2 for  the  given  U2  and  US. 

Compute  the  effective  shear  velocity  from  U\  - V/Ul . 

The  depth  due  to  grain  roughness  o’  = v?/gsa. 

Determine  the  disturbance  coefficient  from  the  following 
procedure: 


(i)  Discharge  and  slope  adjustment: 

(a)  Choose  from  Fig.  8 the  suitable  design  curves  which 
bracket  the  given  q and  S . In  this  example, 

Figs.  8d  and  8e  are  chosen. 


A-l 


(b)  From  Fig.  8d,  a = 0.145  for  q = 50  (ft3/s)/ft, 

F = 0.275,  T = 5°C  and  S = 0.00015;  and  c = 0.044 

o 

for  q = 10  (ft3/s)/ft,  F = 0.275,  T = 5°C  and 
So  = 0.00015.  Then  by  interpolating  on  a loga- 
rithmic graph  of  q versus  a,  we  can  find  a = 0.123 
for  q = 40  (ft3/s)/ft,  F = 0.275,  T = 5°C,  and 
S = 0.00015. 


(c)  Using  the  same  procedure  as  (b) , we  can  find 
a = 0.445  for  q = 40  (ft3/s)/ft,  F = 0.275, 

31  = 5°C,  and  S = 0.0003  from  Fig.  8e. 


(d) 


By  interpolating  on  a logarithmic  graph  of 
a versus  S,  we  can  find  that  a = 0.210  for 
q = 40  (ft3/s)/ft,  F = 0.275,  T = 5°C,  S = 0.0002 
and  = 0.20  mm. 


(ii)  Sediment  grain  size  adjustment 


(a)  Using  Fig.  9,  find  locations  of  a for  ( q = 50 

(ft3/s)/ft,  = 0.2  mm)  and  ( q = 50  (ft3/s)/ft, 

dQn  = 0.22  mm  given)  and  measure  the  vertical  dis- 

01/ 

tance  between  them.  In  this  example,  the  distance 

d » -0.2  cm  (because  c for  0.22  mm  is  below  that 
c 

for  0.20  mm). 

(b)  Measure  the  vertical  distance  between  locations  of 
for  dm  = 0.2  mm  and  discharges  which  bracket  the 

given  discharge  from  Fig.  9.  In  this  example, 

the  distance  d , - 3.72  cm. 
cl 

(c)  Measure  the  distance  between  locations  of  a for 

discharges  which  bracket  the  given  discharge  and 
the  Froude  number  assumed  at  step  1 from  the 
figure  with  slope  closest  to  the  given  slope.  In 
this  example,  the  distance  d - 3.92  cm  from 
Fig.  8d.  c2 

(d)  Compute  the  correction  distance,  d' , = 3.72. 

c 

dc  - (^)OW(^i) 

= (-0.2) (3.92)7(3.72)=  - 0.21  cm. 


A- 2 


(e)  Locate  the  value  of  <?  from  (i)  on  a logarithmic 
graph  with  scale  the  same  as  Fig.  9 and  add 

distance  d'  to  obtain  the  value  of  a = 0.197  for 

a 

q = 40  (ft:l/s)/rt,  r,  = 0.0002,  /■’  = 0.275,  T “ 50C, 

and  d = 0.22  mm. 

m 

Col.  (9):  Find  / from  Figure  10  for  given  q and  T. 


Col.  (10):  Compute  the  value  of  p from 

o 


P 

o 


0-1  )g  iM 
"Y  erj-0.047 


where 


0 = 
O 


(.s-l)g  d5Q 


+ 0.1  s 


Col.  (11):  Compute  the  shear  velocity  from  Eq.  (62),  i.e., 

n ? = fP0ov 

Col.  (12):  Compute  the  depth  due  to  form  roughness  from 

D"  = U”  2/gSo 

Col.  (13):  The  total  depth  equals 

D = D'  + D" 


Col.  (14):  Tne  calculated  discharge  is  then 

q = VD 

If  q > is  different  from  the  given  flow  discharge,  Steps  (1) 

through  (14)  are  repeated  until  the  calculated  discharges 
bracket  the  given  flow  discharge.  The  flow  depth  is  then 
obtained  based  on  linear  interpolation  with  the  discharge. 
In  general,  the  increase  in  the  Froude  number  decreases  the 
calculated  discharge.  After  interpolation,  the  flow  depth 
is  found  to  be  8.55  ft. 


A-3 


2.  Table  1 


Table  1 Summary  of  Field  Data 


Dune 

s 

*50 

d65 

2’ 

D 

L 

H 

(ft/ 

Source 

Date 

rc; 

(fps) 

(ft) 

(ft) 

(ft) 

iooo  ft; 

(mm; 

(mm;  / 

n 

9/8/66 

21.6 

4.48 

10.11 

210 

3.3 

. 15 

.21 

.23  .110 

1.80 

9/13/66 

22.8 

4.51 

10.33 

205 

2.9 

.15 

.21 

.23  .110 

2.28 

9/19/66 

17.8 

4.58 

9.84 

350 

3.3 

.15 

.21 

.23  .117 

1.41 

9/26/66 

17.8 

4.58 

10.00 

335 

3.4 

.15 

.21 

.23  .116 

1.46 

10/3/66 

14.5 

4.72 

9.58 

420 

2.6 

.15 

.21 

.23  .122 

1.62 

10/10/66 

15.6 

5.00 

9.30 

370 

2.8 

. i5 

.21 

.23  .123 

.96 

10/17/66 

9.4 

5.36 

9.35 

450 

2.7 

.15 

.21 

.23  .141 

.76 

10/24/66 

10.0 

5.11 

9.14 

395 

2.5 

.15 

.21 

.23  .133 

.88 

11/7/66 

6.7 

5.12 

9.24 

570 

2.4 

.15 

.21 

.23  .136 

.94 

8/21/67 

22.2 

4.59 

10.13 

250 

3.8 

.15 

.21 

.23  .112 

1.32 

9/25/67 

17.7 

5.19 

9.56 

370 

4.7 

.15 

.21 

.23  .125 

.58 

a> 

> 

10/24/67 

12.3 

5.01 

9.82 

315 

3.8 

.15 

.21 

.23  .124 

.96 

OS 

11/22/67 

5.0 

5.49 

9.11 

170 

1.9 

.15 

.21 

.23  .161 

.80 

•H 

8/6/68 

25.6 

4.23 

11.70 

.15 

.21 

.23  .109 

a 

o 

9/10/68 

18.9 

4.54 

10.15 

229 

3.2 

.15 

.21 

.23  .114 

1.23 

i/) 

t/> 

•H 

10/7/68 

13.3 

5.34 

9.64 

315 

2.8 

.15 

.21 

.23  .129 

1.27 

z 

12/2/68 

2.2 

5.08 

9.67 

2412 

3.0 

.15 

.21 

.23  .137 

1.16 

9/23/69 

19.5 

5.44 

14.04 

224 

3.4 

.15 

.21 

.23  .117 

2.04 

10/14/69 

10.6 

5.66 

12.65 

308 

3.1 

.15 

.21 

.23  .131 

1.28 

11/4/69 

7.2 

5.52 

11.79 

332 

2.7 

.15 

.21 

.23  .134 

1.37 

12/1/69 

2.8 

5.29 

10.08 

1120 

2.3 

.15 

.21 

.23  .144 

1.08 

2/18/52 

3.9 

4.24 

27.90 

.0544 

.19 

.21  .122 

4/27/53 

12.2 

4.42 

33.00 

.0589 

.19 

.21  .117 

Csl 

•r4 

6/13/53 

26.7 

4.00 

30.10 

.0632 

.24 

.27  .102 

P. 

a 

7/20/53  28.9 

3.54 

27.90 

.0603 

.23 

.25  .100 

• H 
(/> 
t/> 

8/10/53  26.1 

3.61 

27.60 

.0661 

.22 

.24  .105 

•r<  f-1 

</>  Q> 
t/>  > 

6/7/54 

19.4 

4.96 

34.20 

.0691 

.23 

.25  .113 

•H  *H 

z a: 

2/24/55 

1.7 

5.11 

34.00 

.0720 

.25 

.27  .120 

to 

Q> 

6/10/65 

17.0 

2.67 

2.76 

50 

2.8 

.64 

.23 

.25  .110 

1.72 

11/29/65 

5.0 

4.98 

3.69 

.66 

.18 

.20  .149 

ctJ 

P 

12/14/66 

2.0 

4.05 

2.48 

.70 

.20 

.22  .138 

u 

12/13/66 

2.0 

4.66 

2.76 

.71 

.20 

.22  .147 

o 

•H 

OS 

6/11/69 

18.0 

3.62 

5.11 

.69 

.28 

.31  .115 

Data  are  from  Corps  of  Engineers (1969)  and  personal  communication  from 
Mr.  Warren  Mellema,  Corps  of  Engineers,  Omaha  District. 

Data  are  from  Jordan  (1965) 

Data  are  from  Culbertson  and  others  (1972). 
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